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Abstract. We present final results from the first phase 
of the EROS search for gravitational microlensing of stars 
in the Magellanic Clouds by unseen deflectors (machos: 
MAssive Compact Halo Objects). The search is sensitive 
to events with time scales between 15 minutes and 200 
days corresponding to deflector masses in the range 10^^ 
to a few M0 . Two events were observed that are compat- 
ible with microlensing by objects of mass « 0.1 Mq. By 
comparing the results with the expected number of events 
for various models of the Galaxy, we conclude that ma- 
chos in the mass range [lO"*", 0.02] Mq make up less than 
20 % (95 % C.L.) of the Halo dark matter. 

Key words: : Galaxy: halo, kinematics and dynamics, 
stellar content - Cosmology : dark matter, gravitational 
lensing 



1. Introduction 

The presence of large quantities of "dark matter" in spi- 
ral galaxies like our own has been inferred from their flat 
rotation curves (e.g. Primack et al, 1988); the dynamic 
mass of the Galaxy is thought to be 3 to 8 times larger 
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than the visible mass up to 50 kpc from the Galactic Cen- 
ter. From primordial nucleosynthesis, we learn that bary- 
onic dark matter can be up to 10 times more abundant 
than visible matter (depending on the value of Hq): all 
Galactic dark matter could thus be baryonic. A possible 
form would be compact objects too light to burn hydrogen 
(m < 0.07 - 0.1 Mq) (Carr, 1990). 

Here we report results from a search for such un- 
seen compact objects in the Galactic Halo performed 
by our collaboration "EROS" (Experience de Recherche 
d'Objets Sombres) at the European Southern Observa- 
tory at La Silla, Chile. Such objects can be detected 
via the gravitational microlensing effect (Paczyhski, 1986) 
that causes an apparent temporary brightening of stars 
outside our Galaxy as the unseen object passes near 
the line of sight. The magnification factor is given by 
A — {u^ + 2)/[u{u'^ + iy^^] where u is the undeflected 
"impact parameter" of the light ray with respect to 
the unseen object in units of the "Einstein Radius", 
Re = {4:GmLx{l - x)l(?fl'^. Here, m is the deflector 
mass, L the observer-source distance and Lx the observer- 
deflector distance. 

The time scale r for the magnification is the time for a 
Halo object to move through a distance equal to its Ein- 
stein radius; its median value is t ~ 90 days ^Jtvl/Mq for 
L= 55 kpc (source star in the LMC) and a standard model 
of the Halo with a velocity dispersion of 245 km.s^^. 

EROS has conducted two observing programs, one us- 
ing a 16 CCD camera mounted on a 40 cm diameter tele- 
scope to search for short time scale microlensing events 
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(r < a few days), and the other using Schmidt photo- 
graphic plates for longer time scales. Results from the 
first 2 years of CCD data were given in (Aubourg et al, 
1995) while results from all 3 years of Schmidt plate data 
were published in (Ansari et al, 1996). In this paper, we 
present final results from the first phase of the EROS pro- 
gram, corresponding to 4 years of CCD observations and 
the Schmidt plate data. 

2. Observations and data reduction 

During 3 annual periods of about 6 months, 290 usable 
photographic plates of 29 x 29 cm^ have been exposed 
at the ESO Im Schmidt telescope, half with a red fil- 
ter and half with a blue filter. The usable field on a 
plate is 5.2° x 5.2°, centered on a = 5h20mn, S = -68°30' 
(eq. 2000). Exposure times were 1 hour in each colour. 
Apart from the very crowded LMC bar region, our star 
detection efficiency abruptly drops at limiting magnitudes 
of about 20.5 in red and 21.5 in blue. The time sampling 
of the plates makes the program sensitive to microlensing 
event durations ranging from one day to a few months. 
Details arc available in (Cavalier, 1994), (Laurent, 1995), 
(Ansari et al, 1996). 

We have also taken during four annual periods more 
than 19,000 images with our CCD setup. The camera 
(Arnaud et al, 1994) consisted of a mosaic of 16 but- 
table 579 X 400 pixels Thomson THX 31157 CCDs. It 
was mounted on a 40 cm reflector (f/10) refurbished by us 
and the Observatoire de Haute-Provence and had a field of 
0.4° X 1.1°. Eleven CCDs were active in 1991-92 during 100 
days and 15 CCDs for the next three seasons of about 230 
days each. The first three years were devoted to the obser- 
vation of one field in the bar of the LMC {a = 5h23.5mn, 
S = — 69°36', eq. 2000), the last year to one field in the 
center of the SMC (a = 0h50mn, 6 = -73°15', eq. 2000). 
Here, we are sensitive to microlensing durations ranging 
from 15 minutes to a few days on stars brighter than about 
19.5 magnitude in V band. Details are available in (Quein- 
nec, 1994), (Aubourg at al, 1995), (Renault, 1996a), (Re- 
nault et al, 1996b). 

The images are processed using a custom designed 
fast photometric reconstruction software to produce light 
curves. Wc then analyse the 250 000 light curves from the 
CCD data and the 6 10^ light curves from Schmidt plates 
in both colours, corresponding to a total of 5 10^ photo- 
metric measurements in crowded fields. 

3. Data analysis 

In order to isolate microlensing events, we searched for 
positive variations occuring simultaneously in the red and 
blue light curves. CCD and Schmidt plate analyses differed 
in details but used the same properties which allow us to 
separate microlensed stars from variable stars. The most 
important one is uniqueness: because of the low predicted 



optical depth (;$ 5.10~^), the probability of a measurable 

microlensing effect to happen twice on the same star in 
a few years time is negligible. We also require equality of 
the base luminosity before and after the variation. In addi- 
tion to these properties, the set of microlensing candidates 
should be representatively distributed throughout our ob- 
served H-R diagram. Neither analysis uses the shape of 
the light curve, thus making us sensitive to events involv- 
ing multiple lenses or sources, and to events affected by 
the finite source size or additional light from unresolved 
stars (blending). 

To obtain detection efficiencies, we processed through 
the same analysis software a sample of observed light 
curves with a random theoretical microlensing shape su- 
perimposed; blending and finite source size effects are 
taken into account. The finite size effect results essen- 
tially in a large loss of efficiency for lensing objects lighter 
than 10~^ Mq; it is negligible for deflectors heavier than 
lO-''-^ Mq (Renault et al, 1996b). The blending correction 
is about four times smaller than that due to the finite size 
of the source for the CCD data. Its impact on the Schmidt 
plates analysis is not important (Ansari et al, 1996). 

No microlensing event was identified in the CCD data. 
Two events compatible with microlensing were identified 
from the Schmidt plate data. In the hypothesis of mi- 
crolensing, they have amplitudes of 1.0 and 1.1 mag. and 
time scales of 23 and 29 days. Those durations correspond 
to deflectors in the mass range [0.01-1 Mq]. The light 
curves were presented in (Aubourg et al, 1993), (Ansari 
et al, 1996). We have taken several spectra and followed 
the candidates with our CCD camera. The first candidate 
is a Bo star (Bcaulicu et al, 1995) and has shown no subse- 
quent variation in plate or CCD observations over a 3 year 
period. The second candidate is an AO star that exhibits 
a 2.8 day periodic variability suggestive of an eclipsing 
system (Ansari et al, 1995). However, the observed sin- 
gle large increases (w 1 mag.) do not correspond to any 
known variability of such stars; we cannot exclude that 
the detected events are due to new types of variable star. 

4. Expected number of events 

The expected number of events can be calculated once a 

model of the Galaxy is chosen. Several models were simu- 
lated using a combination of a disk and a halo. The disk 
models are exponential with a height scale of 0.5 kpc and a 
length scale of 5.0 kpc. We use a surface density Eq of 50 or 
100 MqPC"^ ; the lowest value corresponds to known mat- 
ter in the disk (stars, stellar remnants and gas), whereas 
the largest value is suggested by maximal disk models and 
includes a dark matter component in the disk. 

For the Halo, we have simulated spherical "standard" 
halos and flattened halos (Evans, 1993). Flattened halo 
models are self-consistent and give simultaneously density 
and velocity dispersion. The flattening is characterized by 
the axis ratio q and the shape of the rotation curve by a pa- 
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rameter /3; we assume q = 0.75 and an asymptotically flat 
rotation curve (/? = 0). For spherical and flattened models, 
we use a core radius of 5.6 kpc (Primack et al, 1988) and a 
Galactic Center distance of 7.9 kpc (Merrifield, 1992). We 
scale the halo density distributions such that the rotation 
velocity near the Sun be 200 km.s"^ (Merrifield, 1992). We 
have also simulated our "reference model" used in (Ansari 
et al, 1996). It is a "standard" spherical model with a core 
radius of 7.8 kpc and a velocity dispersion of 245 km.s~^. 
The mass of its Halo is normalised to 4.10^^ Mq within 
50 kpc. Table 1 gives characteristics for five models: the 
local Disk surface density, the local Halo density and the 
circular velocity at 50 kpc. Other models are described in 
(Renault et al, 1996b). 
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Table 1. Characteristics of the five models; model 5 is the 
reference model described in (Ansari et al, 1996). Eq is the 
local Disk surface density and p© is the local Halo density. 
The circular velocities Vdrc are those generated by the mass 
in the Halo or by the total mass (disk+halo) at a radius equal 
to the LMC distance. 

The expected number of events due to Halo objects 
as a function of their mass is plotted in figure 1, where 
it is assumed that all machos have the same mass m. We 
have combined Schmidt plates and CCD programs. Two 
comparable maxima are seen where observations are most 
efflcient: at m = 3 10~^ M© with the CCD data and 
m = 10~^ Mq with the Schmidt plates. The expected 
numbers of events for the five models approximately scale 
as Vj^c (table 1) and therefore reflect the halo masses. 

The expected number of events due to lensing by 
Galactic stars is less than 0.15; for lensing by stars in the 
Magellanic Clouds, it is less than 0.6. 

5. Constraints on the Halo 

From the expected and observed number of events, we ob- 
tain an upper limit on the fraction of the halo along the 
line of sight that is composed of machos. The statistical 
method is the same as in (Ansari et al, 1996). Figure 2 
presents 95 % CL limits for the reference model, all the 
EROS data and two microlensing candidates. This limit 
is shown as a function of the assumed deflector mass. For 
the CCD data, we also show the effect of taking into ac- 
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Fig. 1. Expected number of microlensing events in the EROS 
programs (CCD -LMC and SMC sources- and Schmidt plates) 
assuming that all Galactic dark matter is in the form of machos 
of the same mass. The five curves refer to the Galajcy models 

of table 1. The full curve corresponds to the reference model, 
dashed curves to fiattened models and dotted curves to the 
spherical models. 



count blending and finite size effects, important only be- 
low 10~^ M©; these effects supress all sensitivity below 
10~^ Mq. We also show limits assuming or 1 microlens- 
ing events, that are identical for objects with masses below 
10"'^ Mq. For two microlensing events, we can not give 
significant constraints beyond 1 Mq. 

Figure 3 presents the limits obtained for different mod- 
els of the Galaxy and all the EROS data. The hypothesis 
that the two observed events arc true microlensing events 
is adopted in order to obtain a conservative upper limit. 

The used disk length scale (DLS) is likely to be high 
(Sackett, 1996) but models with low disk surface density 
and a DLS of 2.5 kpc are more constraining than our mod- 
els with a heavy disk and a DLS of 5 kpc. Moreover, heavy 
disk and low DLS (<3.2 kpc) are not compatible with the 
observed rotation curve. 

We can divide the studied mass range [10~^-1 Mq] in 
three distinct parts. 

1) No limits can be inferred for objects lighter than 
10"'' Mq; the very existence of such objects is doubtful. 

2) For an intermediate model, we observe that the frac- 
tion of the Halo in the form of objects with masses be- 
tween lO"'' and 0.02 Mq is below 20 % (between 5 10"'' 
and 0.002 Mq, it is lower than 10 %). As the limit is rather 
mass independent, it is valid for any mass function in this 
interval. This limit is similar to that of the MACHO col- 
laboration (Alcock et al, 1996a) but extends an order of 
magnitude lower in mass. We note that the limits from 
the two groups are derived from almost independent sets 
of stars observed at different times. Thus they could be 
combined to obtain an even lower limit. 

3) The limit for objects with masses between 0.02 and 
0.5 Mq becomes less and less stringent at higher masses. 
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Fig. 2. Exclusion diagram at 95 % CL for the reference model 
with all EROS data, assuming all deflectors to have the same 
mass. For the CCD program, wc show the influence of blending 
and finite size effects (the dotted line on the left is the limit 
without those effects). Limits are shown for (dashed line), 
1 (mixed line) or 2 (full line) candidates assumed to be ac- 
tual microlensing. The cross is centered on the area allowed at 
95 % CL by the MACHO program (Alcock et al, 1996b) as- 
suming 6 microlensing events and a standard spherical model 
(model S, very similar to our reference model). We also indi- 
cate the MACHO exclusion contour obtained by combining all 
their results with t < 20 days (Alcock et al, 1996a). 

No robust limits can be inferred for objects heavier than 
0.2 because, with our present sensitivity, some models 
are compatible with a Halo entirely made of such objects. 
Our results do not contradict the positive signal of the 
MACHO collaboration (Alcock et al, 1996b) as is shown 
in figure 2. 

6. Conclusions 

The main result of the first phase of the EROS program 
is that objects with masses between 10"'^ and 0.02 M© 
do not contribute significantly to the Halo dark matter. 
The results shown in figure 2 indicate that there remains a 
great uncertainty in the proportion of the Halo comprised 
of heavier objects. Further data from the MACHO col- 
laboration and the recently upgraded OGLE and EROS 
programs should clarify the issue in the coming years. 
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Fig. 3. Exclusion diagram at 95 % CL assuming the same mass 
for all deflectors, for 5 Calaxy models of table 1. The full curve 
corresponds to the reference model, dashed curves to flattened 
models and dotted curves to the spherical models. We have 
used all EROS data and assumed two observed microlensing 
events with r = 23 and 29 days. 
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